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Section  1 
INTRODUCTION 


Measurements  from  continuously  recording,  in  situ, 
sound  velocity,  salinity,  temperature,  depth  probes  TSvsfDs) 
at  262  stations  are  used  to  describe  the  distribution  of 
mixed  layer  depths  (MLDs)  in  the  northern  North  Atlantic. 
Because  the  term  mixed  layer  refers  to  a  sea  surface  layer 
with  nearly,  but  not  completely,  homogeneous  temperature 
and  salinity  properties,  several  definitions  of  MLD  are 
possible.  Four  definitions  based  on  four  parameters  are 
used  here.  The  four  parameters  are  temperature,  sound 
speed,  density  and  stability. 

The  purpose  of  this  work  is  twofold.  The  first 
purpose  is  to  determine  if  there  are  any  significant  differ¬ 
ences  in  the  distribution  of  MLDs  among  various  seasonal  and 
spatial  subsets  of  stations  in  the  area  of  interest.  Eleven 
such  subsets  (called  "classes")  have  been  identified.  The 
second  purpose  is  to  determine  if  there  are  any  significant 
differences  in  the  distribution  of  MLDs  due  to  the  four 
different  definitions  of  MLD. 

The  vertical  profiles  from  the  SVSTD  measurements 
are  smoothed  to  reduce  the  effect  of  noise  and  the  effect  of 
f i nestructure  with  vertical  wave  lengths  smaller  than  five 
meters  (5m).  The  details  of  this  smoothing  process  are 
presented  in  section  2.1.  The  four  parameters  used  in  the 
definitions  of  MLD  are  calculated  for  each  vertical  level 
using  the  algorithms  detailed  in  section  2.2.  In  that 
section,  also,  the  distinction  between  depth  and  pressure, 
which  Is  maintained  throughout  the  report,  is  elucidated. 

In  section  3,  the  definitions  of  MLD  based  on  each 
of  the  four  parameters  are  stated  and  reference  is  made  to 
the  computer  routines  written  to  implement  the  definitions. 
Listings  of  the  computer  routines  are  provided  in  Appendix 
A.  The  eleven  classes  are  then  defined.  Finally  presented 
In  section  3  are  the  results:  the  distributions  of  MLD  in 
the  form  of  histograms  for  all  eleven  classes  for  each  of 
the  four  definitions. 

Each  distribution  is  characterized  in  section  4  by 
its  mean  and  its  variance.  Confidence  intervals  for  these 
parameters  are  presented  based  on  the  assumption  that  the 
MLDs  represent  samples  from  a  Gaussian  distribution.  The 
degree  to  which  the  assumption  of  normality  is  justified  is 


tested  with  a  chi-square  goodness-of-f it  test.  Section  4 
also  identifies  the  significant  variations  of  the  MLD 
distribution  characteristics  with  class  and  def  1  n  1 1  i  <5n 
through  the  use  of  non-par  am etric  statistical  methods. 

The  major  results  of  this  investigation  are  sum¬ 
marized  in  section  5.  Items  requirlnq  further  investigation 
are  identified  in  section  6. 


Section  2 

PROFILE  PREPARATIONS 


The  SVSTO  stations  are  preconditioned  by  smoothing 
vertical  profiles  of  temperature  and  salinity  and  then 
calculating  values  of  sound  speed,  density  (specific  gravity 
anomaly)  and  stability  { Brunt-VSi s‘41*4  frequency). 


2.1  FILTER 

The  data  are  supplied  with  a  resolution  of  1.0047 
dbar,  but  noise  in  the  data,  including  quantizing  noise 
introduced  by  rounding  both  temperature  and  salinity  values 
to  the  nearest  hundredth,  has  a  large  effect  on  the  vertical 
gradients  calculated  by  taking  differences  over  such  small 
scales.  Consequently,  filtering  was  performed  to  reduce 
the  contribution  to  the  stratification  by  components  with 
vertical  scales  less  than  5  dbar.  Because  many  of  the 
mixed  layer  depth  definitions  require  a  comparison  to  a 
value  at  5m  and  because  the  shallowest  values  at  many 
stations  is  5m,  the  filtering  scheme  must  retain  end 
values . 


In  an  ideal  situation,  the  first  observation 
is  at  0  dbar  so  an  eleven  weight  filter  would  return  a 
value  at  nominally  5  dbar.  For  this  reason  only  filters 
with  eleven  weights  or  less  were  considered.  Using  a 
computer  program  that  will  design  finite  impulse  response, 
linear  phase,  digital  filters  using  the  Remez  Exhange 
Algorithm  (McClellan,  1973),  an  eleven  point  filter  was 
designed  that  would  pass  components  with  scales  larger  than 
10  dbar  and  would  strongly  attenuate  (i.e.,  stop)  components 
with  scales  less  than  5  dbar.  In  the  filter  design  there  is 
an  inverse  relationship  between  the  degree  of  attenuation  in 
the  stop  band  and  the  degree  of  uniformity  of  gain  in  the 
pass  band.  After  trying  several  cases  with  a  nine  weight 
filter  It  was  determined  that  nonuniformity  in  the  gain  of 
the  pass  band  did  not  affect  the  large  scale  features  of  the 
stratification  noticeably.  That  is,  the  location  and  size 
of  large  scale  variations  did  not  show  shifts  of  a  decibar 
or  more.  The  degree  of  attenuation  in  the  stop  band  did 
have  a  noticeable  effect  on  the  gradient  dependent  param¬ 
eters  such  as  Brunt-VSi  s SCI  If  frequency.  Therefore,  in  the 
filter  design,  emphasis  was  placed  on  attenuating  the  short 
wave  lengths. 
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The  smoothed  profiles  produced  by  the  nine  weight 
filters  were  compared  to  those  of  an  eleven  weight  filter. 
Use  of  the  eleven  weight  filters  did  not  result  in  a  notice¬ 
ably  smoother  Brunt- V'4i  s^lli  profile;  therefore,  the  extra 
calculations  necessary  for  the  eleven  weight  filter  were 
deemed  unwarranted.  The  nine  weight  filter  was  selected 
for  smoothing  the  profiles  of  temperature  and  salinity.  The 
weights  of  the  chosen  nine  point  filter  are  given  in  table 
2.1,  the  impulse  r,esponse  is  given  in  figure  2.1a  and  the 
wave  number  response  is  given  in  figure  2.1b. 

In  order  that  the  filter  may  be  applied,  the 
input  data  must  consist  of  a  uniform  pressure  series. 
Consequently  linear  interpolation  was  performed  to  generate 
values  of  temperature  and  salinity  every  1.0047  dbar. 

In  order  to  retain  the  end  values  of  the  pressure 
series,  filters  with  odd  numbers  of  weights  less  than  nine 
must  be  applied.  The  nine  point  filter  gives  a  value 
appropriate  at  the  fifth  pressure  in  the  series,  a  seven 
point  filter  provides  a  value  appropriate  at  the  fourth,  a 
five  point  filter  at  the  third,  and  a  three  point  filter  at 
the  second.  The  value  at  the  first  pressure  is  retained  at 
its  original  magnitude.  This  scheme  is  applied  to  the 
beginning  of  the  pressure  series  and  the  weights  for  these 
filters  are  given  in  table  2.2. 

Profiles  of  temperature,  sound  speed,  density  and 
Brunt-Vai sal  a  frequency  (stability)  at  a  sample  station  are 
shown  for  both  original  data  (figures  2.2  and  2.3)  and  for 
smoothed  data  (figures  2.4  and  2.5). 


2.2  ALGORITHMS 

Mixed  layer  depths  are  to  be  defined  in  terms 
of  several  parameters.  All  these  parameters  can  be  cal¬ 
culated  from  pressure  series  of  temperature  and  salinity. 
This  section  presents  the  formulas  for  the  calculation  of 
the  parameters  of  interest:  depth  (D),  sound  velocity  (C), 
sigma-T  (a^)  and  Brunt- Vai sal  a  frequency  ( N ) . 


2.2.1  Pressure-Depth 


Pressure  (P)  measured  by  the  STD  was  converted  to 
depth  (Z)  by  N00  personnel  (R.  Rushton,  personal  communica¬ 
tion)  using  the  relation 


Z 


P 

1.0047 


(2-1) 
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TABLE  2.1 


NINE  POINT  FILTER  WEIGHTS 


POSITION 

WEIGHT 

1 

.00954717 

2 

.06623742 

3 

.12782102 

4 

.18947470 

5 

.21383940 

6 

.18947470 

7 

.12782102 

8 

.06623742 

9 

.00954717 

(NJ  (O  <  IT)  VO  N 


TA8LE  2.2 


WEIGHTS  FOR  THE  SEVEN,  FIVE,  THREE 
AND  ONE  POINT  FILTERS 
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PRESSURE  (DB) 


PRESSURE  (DB) 


SOUND  SPEED  (M/SEC) 

1490  1500 


where  Z  is  depth  (m)  and  P  is  pressure  (dbar).  Pressure  is 
desirable  for  many  calculations  (sound  velocity  and  in  situ 
specific  volume)  so  the  vertical  coordinate  was  transTormed 
back  to  pressure  using  the  inverse  of  equation  (2-1)  before 
filtering.  Oepth  at  the  nth  observation  can  then  be  calcu¬ 
lated  from  pressure  and  j_n  situ  specific  volume  by  perform¬ 
ing  the  summation  based  upon  the  hydrostatic  relation 

n 

z  -  104  E  ®i(pi"pl-l )  ^Wl)  (2-2) 

i*l 

where 

Pi  =  pressure  at  level  i  (dbar),  P0  *  0 
ai  =  III  s j tu  specific  volume  at  level  i, 

‘  °o  •  «1 

■j  =  gravitational  acceleration  at  level 

i  (m  sec"*) 

The  value  of  gi  can  be  calculated  as  by  Saunders  and  Fofo- 
noff  (1976) 

g.  =  9. 780318  (l+5. 3024x1 0_3sln24>-5.9xl0~6sin224)) 

-6  U-3) 

+  1.113  x  10"° 

where  d>  is  latitude. 

The  value  of  is  determined  from  Ti  and  Si  using  equation 
(2-15). 


2.2.2  Sound  Speed 


Sound  speed  depends  upon  pressure,  temperature 
and  salinity.  To  be  most  consistent  with  NAVOCEANO  cal¬ 
culations  of  sound  speed  the  Wilson  (1960)  formula  is  used 
in  this  study.  Wilson  uses  a  pressure  unit  of  kg/cm^. 
Pressure  in  the  ocean  is  expressed  in  terms  of  decibars 
above  atmospheric  pressure.  The  two  units  must  be  related 
to  use  the  formulas. 

If  PO  is  the  oceanographic  pressure  in  dbar,  and 
P  the  total  pressure  in  kg/cm‘,  then 

p  =  (2-4) 
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This  pressure,  plus  temperature. (already  In  °C)  and 
salinity  (already  in  o/oo )  ,  are  used  In  the  formula 
of  Wi 1  son 


C  -  SV(T,S,P)  »  1449.14  +  Vy  +  Vp  +  Vs  +  VSyp  (2-5) 

where 

Vy  «  4.5721  T  -  4.4532xl0“2  T2  -  2.6045xl0"4  T3  +  7.9851x10-6  y4 
Vp  «  1 .60272x10-1  P  +  1.0268x10-5  p2  +  3.521x10-9  p3  -  3.3603x10-12  P4 
VS  -  1.39799(S-35)  +  1 .69202xl0-3(S-35)2 
and 

V$Tp  *  (S-35)  (-1. 1244xl0-2  T+7 . 7711 xlO"7  T2+7. 701 6x10*5  P 
-1. 2943xl0-7  P2+3. 158x10-8  py+1 .5790x10-9  py2) 

+P(-1.8607xl0-4  T+7. 4812x10-6  T2+4. 5283x10-8  y3) 

+P(-2. 5294x10-7  T+l .8563x10-9  y2) 

+p3(-l. 9646x10-10  T) 


2.2.3  Siqma-T 

The  specific  gravity  anomaly,  ot,  is  calculated 
using  the  formula  of  Knudsen  (1901)  with  a  slight  modifica¬ 
tion  In  the  salinity  (S)  -  chlorinity  (Cl)  relationship. 

Whereas  Knudson  defined  S  «  0.03  +  1.805  Cl,  in  this  work 
the  following  definition  is  used: 

S  *  1 .80655  Cl .  (2-6) 

This  ensures  that  both  salinity  and  chlorinity  are  con¬ 
servative  properties.  Equation  (2-6)  is  then  used  to  define 
a  chlorinity  value  to  insert  in  the  Knudson  formula 

o y(T,S)  »  Zy  +  (oQ+.  1324)  (l-AT+BT(o0-.  1324))  (2-7a) 

where 

oo«.-069  +  1.4708  x  Cl  -  .001570  x  Cl2  +  .398  x  10-4  x  Cl3  (2-7b) 
£y  -  -(T-3.98)2(T+2.83)/ [(503. 57)  (T+67. 26)] 

At  -  T  (4. 7867- . 0981 85T+.  1 0843x1 0"2T2)  x  10‘3 
Bt  «  T  (l 8. 03-  .  81 64T+.  01 667T2)  x  10"6 
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where 

12  =  depth  of  deeper  observation 

Zi  =  depth  of  shallower  observation 

=  density  of  deeper  observation  at  a 
reference  pressure 

Pj  *  density  of  the  shallower  observation 
at  the  same  reference  pressure. 

The  density  in  equations  (2-8)  and  (2-9)  is  computed  from  an 
equation  by  Ekman  (1908)  for  in  situ  specific  volume  where 
a  reference  pressure,  PR,  midTway  ETitween  the  two  observa¬ 
tions  is  used.  The  density  equation  depends  upon  the 
reference  pressure,  the  salinity  of  the  observation  and  the 
potential  temperature  (®pr)  of  the  observation  at  the 
reference  pressure.  6  can  be  computed  from  T,  the 
pressure  difference  (PR-P)  and  the  adiabatic  lapse  rate 
(y  )  .  The  value  of  Y  can  be  computed  from  a  formula  by 
Bryden  (i97i)  which  depends  on  temperature,  salinity  and 
pressure.  The  pressure  used  for  the  calculation  of  6  pr  is 
a  pressure  midway  between  the  observed  pressure  and  the 
ref or  20' 2  0  e *  s j r e  . 
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The  algorithm  for  the  computation  of  N  from 
observations  of  Pi,  Ti,  Sj,  and  P2,  T2,  S  2  therefore 

consists  of  the  following  steps. 

PR  *  (Pt  +  P2)  (2-10) 

PPj  «  (Pj  +  PR)  for  i  *  1,2  (2-11) 

From  Bryden  (1973), 

Y1  «  10-3  SEE  A i kc  P P i j  (  S i  -  35  ) k  t/  (2-12) 

3  lc  £  J  * 


for  i  =  1,2 


where 


A000 

s 

.35803 

X 

10-1 

A1 02 

- 

.87330 

X 

10-8 

A001 

s 

-85258 

X 

lO-2 

A 103 

3 

-.54481 

X 

0 

1 

0 

A002 

s 

-.68360 

X 

10*4 

A1  10 

S 

-.11351 

X 

10-6 

a003 

s 

.66228 

X 

10-6 

A1  11 

3 

.27759 

X 

00 

1 

0 

A010 

= 

.18932 

X 

10-2 

a2C0 

3 

-.46206 

X 

10-9 

A01 1 

s 

-.4 2393 

X 

10-4 

A  2  0 1 

3 

.  18676 

X 

10-10 

A100 

3 

.18741 

X 

1 

O 

a202 

S 

-.21687 

X 

10-12 

Then, 


0phi  *  Ti  +  (re-pi)  n 


for  i  »  1,2 


-i  -  .(!*,! t,»J-  to  1-1,2  (2_14) 

From  Ekman  (1908), 

“(^RI.Si.PR)  -  cto1-[PR(aoiXio-9)(A-Bi+Ci-Di-Fi+Gi)l  (2-15) 


( 2-13) 

(2-14) 


where 


4886 

1  +.  183x1  0"4Pr) 


3i 

=  227  +  28 

.  336 

PRi 

-551  52r. 

+ 

.004  9 

ci 

=  PR  10‘4( 

105. 

5  +  9. 

50  e  pRi - . 1 58 

6  PRi) 

Di 

=  1.5  PR2 

9  PRi 

10" 

8 

Fi 

-  /  a  -28  ) 

■  (  ?o  1 

1  (l  47 

.3-2 

.72  epRi+0. 

04 

-PR 

10'4  (32.4- 

87 

0  PRi 

1  *  (4.S-0-!  epRi-PR  TO-4  (1.8-0.6  epR  )) 

a  _  ,  _ ] _ 

( 1+aT  (e PR i  ,Si )  10  3) 


Both  and  aj  (9  p  R  i  ,  s -j  )  ,  are  defined  according  to  the 

relations  of  Knudsen  (1901)  given  in  equations  ( 2 - 7  a )  and 
(2-7b). 

Finally,  the  depth  difference  is  approximated  by 
the  pressure  difference 


(Z2  -  l])  *  (P 2  -  pl)  (2-16) 

The  values  of  p ] ,  P2  and  (Z2  -  Z i )  are  then  substituted 
into  equation  (2-9)  to  produce  a  value  of  N. 


2.3  ERROR  ESTIMATES 

Estimates  of  errors  in  the  four  parameters  used  to 
define  mixed  layer  depths  are  necessary  if  the  values  are  to 
be  used  properly.  The  uncertainty  in  all  parameters  depends 
upon  the  uncertainties  in  the  basic  parameters  pressure, 
temperature  and  salinity.  The  original  data  are  reported 
every  1.0047  dbar,  to  a  resolution  of  .01°C  and  .01  °/oo. 
These  values  can  be  taken  as  the  precision  of  the  tempera¬ 
ture  and  salinity,  respectively.  The  smoothed  profiles 
consist  of  values  still  every  1.0047  dbar  but  averaged  (with 
uneven  weights)  over  nine  points.  The  major  contribution  to 
the  average  comes  from  the  central  three  weights  so  that  for 
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the  purpose  of  approximating  the  precision  error  in  the 
smoothed  temperature  and  salinity  an  average  over  3  points 
can  be  assumed.  Then  the  uncertainty  in  temperature  is 
•01°C/  3  *  .006°C  and  the  uncertainty  in  salinity  is 
.01  °/oo/  3  *  .006  °/oo.  Such  uncertainties  in  temperature 
and  salinity  affect  the  calculation  of  sound  speed  by  less 
than  .05m  sec"^  over  the  range  of  temperatures ,  salinities 
and  pressures  encountered  in  this  work.  The  effect  of 
errors  in  temperature  and  salinity  on  o  y  is  .01  v  y  units. 
Therefore,  changes  in  cr  y  can  be  estimated  with  an  error  of 
under  .015  Qj  units.  The  effect  on  Brunt-Vaisal’a  frequency 
(N)  requires  a  more  complicated  calculation  and  is  found  to 
be  a  function  of  the  value  of  N.  At  values  of  1  cycle 
per  hour  (cph),  the  error  is  about  lcph.  At  2cph  the 
error  is  .5cph  and  at  4cph  the  error  is  .25cph. 
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Section  3 

DISTRIBUTIONS  OF  MLD 


The  smoothed  profiles  of  temperature  (T),  sound 
speed  (C),  specific  volume  anomaly  (cy),  and  Brunt-Vai sal  a 
frequency  (N)  are  used  with  four  definitions  of  mixed  layer 
depths  to  produce  MLDy,  MLDq,  MLDp.  and  MLDn,  respec¬ 
tively.  Each  of  the  262  stations  in  the  northern  North 
Atlantic  is  assigned  to  several  of  eleven  total  classes 
based  on  season  (winter,  spring,  summer  and  fall)  and  sub- 
regions  (termed  areas  I,  II  and  III).  The  distribution  of 
each  of  the  four  kinds  of  MLD  for  each  of  the  eleven  classes 
are  presented  in  this  section. 


3.1  FOUR  DEFINITIONS  OF  MIXED  LAYER  DEPTH 

In  all  these  definitions  the  first  5m  (5.0235 
dbar)  of  data  are  excluded. 


3.1.1  Temperature  Mixed  Layer  Depth  (MLDy) 


According  to  the  temperature  criterion  the  mixed 
layer  is  defined  as  the  near  surface  depth  at  which  the 
absolute  value  of  the  temperature  gradient  exceeds  ,005°C 
m-1  and  at  which  the  temperature  itself  differs  by  at  least 
0.2°C  from  that  at  5m.  Applying  this  criterion  to  the 
pressure  series  .005°C  m""1  was  taken  to  be  equivalent  to 
.005°C  per  1.0047  dbar.  The  computer  subroutine  that 
implements  this  criterion  is  provided  in  appendix  A,  section 
1  . 


3.1.2  Sound  Speed  Mixed  Layer  Depth  (MLD^) 


Accord  ':nq  to  the  sound  speed  criterion  the  mixed 
layer  is  defined  as  the  near  surface  depth  at  which  the 
first  maximum  in  sound  velocity  occurs.  During  the  imple¬ 
mentation  of  this  criterion  a  maximum  was  not  considered 
unless  it  was  greater  than  deeper  values  of  sound  speed  by 
.5m  sec"^  which  is  more  than  the  expected  error.  Before 
this  test  the  sound  speeds  are  rounded  to  the  nearest  .lm 
sec-'.  The  computer  subroutine  that  implements  this  cri¬ 
terion  is  provided  in  appendix  A,  section  2. 
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According  to  the  density  criterion  the  mixed 
layer  is  defined  as  the  near  surface  depth  at  which  the 
first  maximum  positive  density  gradient  occurs  and  at  which 
the  density  itself  exceeds  that  at  5m  by  at  least  .02 
sigma-T  units.  Ouring  imp  1 ementat i on ,  density  data  within 
.02  sigma-T  units  are  skipped  first.  Once  large  enough 
density  values  are  found  a  maximum  difference  is  con¬ 
sidered  if  it  is  greater  than  deeper  differences  by  .01 
sigma-T  units  or  more.  Though  the  exp.ected  noise  in  the 
difference  is  .015  sigma-T  units,  the  described  implementa¬ 
tion  gives  results  that  compare  favorably  to  choices  made  by 
visually  inspecting  the  density  profile.  The  computer 
subroutine  that  implements  the  criterion  is  provided  in 
appendix  A,  section  3. 


3.1.4  Stability  Mixed  Layer  Depth  ( ML ) 


According  to  the  stability  criterion  the  mixed 
layer  is  defined  as  the  near  surface  depth  at  which  the 
first  maximum  exceeding  2  cycles  per  hour  (cph)  in  the 
Brunt-Vai sal  a  frequency  (N)  occurs.  During  implementation  N 
values  less  than  2cph  are  not  considered  until  a  value 
greater  than  2cph  is  found.  The  maximum  is  then  defined 
if  a  deeper  value  of  N  is  less  than  the  maximum  by  more 
than  .5cph  which  is  the  expected  error  at  2  cph.  The 
computer  subroutine  that  implements  this  criterion  is 
provided  in  appendix  A,  section  4. 


3.2  ELEVEN  CLASSES  OF  STATIONS 

Data  provided  for  this  study  consist  of  stations 
that  span  several  seasons  in  three  subregions  termed  I 
(Rockall  Trough),  II  (Iceland  Flank)  and  III  ( I  cel and-Faeroe 
Frontal  Zone).  Figure  3.1  shows  a  map  locating  these  three 
areas.  In  area  I  data  is  available  in  summer  and  fall,  in 
area  II  data  is  available  in  summer  and  fall,  and  in  area 
III  there  is  data  in  spring,  summer  and  fall.  Each  season 
for  each  area  constitute  each  of  the  first  seven  classes. 
The  cruises  and  stations  making  up  these  basic  classes  are 
summarized  in  table  3.1.  The  remaining  classes  consist  of 
all  seasons  for  each  area  (three  classes)  and  all  seasons 
for  all  areas  (one  class). 
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TABLE  3.1 


STATIONS  PROVIDED  FOR  STUDY 


AREA 

SEASON 

CRUISE 

STATION  NUMBERS 

I 

Summer 

932014 

343406 

343405 

1-3,  20-29 

20-44,  46,  47 

39-44 

I 

Fa  1 1 

933005 

57-66 

II 

Summer 

932014 

343406 

343405 

4-19 

1-19,  49-55 

1-24,  26-38,  55-58 

I  I 

Fal  1 

933005 

41-43,  48-54 

III 

Spr i ng 

343525 

1-2,  17-20,  22, 

24-36,  38-47,  50-54 

III 

Summer 

343405 

343622 

45-54 

1-56 

I  I  I 

Fal  1 

343606 

2-4,  6,  17-19,  22, 

24,  29,  31,  36,  37 

j 
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The  computer  tapes  provided  by  NAVOCEANO  also 
contained  additional  stations,  duplicate  stations  and 
several  additional  casts  of  each  of  two  yo-yoed  stations. 
All  these  data  were  iqnored  in  compiling  the  statistics 
presented  here. 


3.3  HISTOGRAMS  OF  MLDs 

The  four  mixed  layer  depths  found  for  each  station 
were  assembled  in  the  eleven  classes  just  described  to  give 
44  distributions  of  MLD .  These  distributions  are  presented 
here  as  histograms  of  percent  frequency  of  occurrence  per  5 
meter  bin  for  the  upper  two  hundred  meters  in  figures  3.1  to 
3.44.  The  first  bin,  includinq  occurrences  in  the  depth 
interval  0  to  4.9999m,  can  have  no  members  because  all 
four  definitions  of  MLD  excluded  the  first  5m  of  data. 
For  these  histograms  depth  was  not  determined  from  pressure 
by  the  numerical  integration  of  equation  2-2.  Instead  depth 
values  in  meters  were  considered  identical  to  pressure 
values  in  decibars  rounded  to  the  greatest  integer.  This 
simplification  saves  many  computations  and  is  only  likely  to 
be  in  error  by  a  fraction  of  a  meter.  All  mixed  layer 
depths  greater  than  200m  are  assigned  to  the  last  bin  in  the 
histograms.  MLDj  is  presented  in  figures  3.2  to  3.12, 
MLDc  in  figures  3.13  to  3.23,  MLDg  in  3.24  to  3.34  and 
MLDn  in  figures  3.35  to  3.45. 

The  salient  feature  of  the  cumulative  histograms 
(figures  3.12,  3.23,  3.34  and  3.45)  are  that  sound  speed  and 
stability  criteria  give  distributions  of  MLD  that  peak  in 
the  shallowest  allowable  bin  while  both  temperature  and 
density  criteria  give  distributions  that  peak  at  depth  (20- 
25m  MLDj  and  30-35m  for  MLDn).  The  cause  for  this  may 
be  that  the  latter  two  have  a  stipulation  in  their  defini¬ 
tions  that  the  value  at  the  base  of  the  mixed  layer  differ 
from  the  5m  value  by  a  minimal  amount.  The  first  two 
definitions  do  not  specify  a  minimum  difference. 

The  summer  histograms  (figures  3.2,  3.3,  3.4, 
3.13,  3.14,  3.15,  3.24,  3.25,  3.26,  3.35,  3.36,  3.37)  show 
that  area  II  has  the  deepest  MLDs  and  the  greatest  variance 
for  all  criteria.  While  both  area  I  and  area  III  have 
shallower  MLDs,  area  III  shows  the  smallest  variance  for 
MLD-j-  and  MLDq.  For  MLDc  and  MLDn  area  I  is  more  stronqly 
peaked  in  the  shallowest  bin  (5-10m).  Fall  distributions 
are  deeper  than  summer  for  all  criteria  but  the  deepest  MLDs 
are  found  for  a  small  percentage  of  the  spring  observations. 


Map  showing  the  locations  of  the  three  areas 
In  which  the  SVSTD  stations  are  grouped. 


Figure  3.3  Histogram  of  MLDT  for  Summer,  Area  II  (77  Stations) 


Figure  3.4  Histogram  of  MLDT  for  Summer,  Area  III  (63  Stations) 
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Figure  3.6  Histogram  of  MLDT  for  Fall,  Area  II  (9  Stations) 
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Figure  3.9  Histogram  of  MLDT  for  All  Seasons,  Area  I  (54  Stations) 


Igure  3.11  Histogram  of  MLDT  for  All  Seasons,  Area  III  (106  Stations] 


Histogram  of  MLD,  for  All  Seasons,  All  Areas  (246  Static 
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Histogram  of  MLDr  for  Summer,  Area  I  (44  Stations) 


Histogram  of  MLDr  for  Summer,  Area  II  (77  Stations) 
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Figure  3.16  Histogram  of  HLOr  for  Fall,  Area  I  (9  Stations) 


Figure  3.17  Histogram  of  MLD-  for  Fall,  Area  II  (9  Stations) 


Figure  3.18  Histogram  of  #IDC  for  Fall,  Area  III  (10  Stations) 


re  3.19  Histogram  of  MLDC  for  Spring,  Area  III  (31  Stations) 


gure  3.21  Histogram  of  MLD-  for  All  Seasons,  Area  II  (86  Stations) 


Figure  3.22  Hist 
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Igure  3.26  Histogram  of  MLDg  for  Summer,  Area  III  (63  Stations) 


Igure  3.27  Histogram  of  MLDn  for  Fall,  Area  I  (9  Stations) 


gure  3.28  Histogram  of  MLD0  for  Fall,  Area  II  (9  Stations) 


gure  3.29  Histogram  of  MLDp  for  Fall,  Area  III  (11  Stations) 


gure  3.31  Histogram  of  MLDD  for  All  Seasons,  Area  I  (54  Stations) 


gure  3.32  Histogram  of  MLDd  for  All  Seasons,  Area  II  (86  Stations) 


Figure  3.34  Histogram  of  MLDD  for  All  Seasons,  All  Areas  (246  Station 


Figure  3.35  Histogram  of  MLDN  for  Summer,  Area  I  (45  Stations) 


Mgure  3.36  Histogram  of  MLDN  for  Summer,  Area  II  (77  Stations) 
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Figure  3.38  Histogram  of  MLDM  for  Fall,  Area  I  (9  Stations) 


Figure  3.39  Histogram  of  MLDW  for  Fall,  Area  II  (9  Stations) 


32  Stations) 


gure  3.43  Histogram  of  MLDN  for  All  Seas 


1 


65 


gure  3.45  Histogram  of  MLDN  for  All  Seasons,  All  Areas  (246  Stat 


Section  4 

SIGNIFICANT  VARIATIONS 


The  information  in  the  large  number  of  histograms 
presented  in  section  3  is  summarized  here  and  analyzed  by 
means  of  both  parametric  and  non-par ametri c  statistical 
methods  to  indicate  the  most  significant  relationships. 


4.1  MEANS  AND  STANDARD  DEVIATIONS  OF  MLDS 

The  histograms  presented  in  the  previous  section 
can  be  characterized  by  their  means  and  standard  deviations. 
These  values  are  presented  in  table  4.1  along  with  estimates 
of  their  confidence  intervals  based  upon  Student-t  and 
chi-square  distributions.  Both  these  computations  assume 
that  the  underlying  distributions  are  Gaussian,  an  assump¬ 
tion  that  is  not  ideal,  as  will  be  shown.  However,  the 
intervals  are  useful  as  first  order  estimates  of  statistical 
significance. 

The  entries  in  table  4.1  are  summarized  in  figures 
4.1  to  4.4  which,  for  each  definition  show  how  the  mean  MLD 
varies  with  season  and  area.  Similar  plots  are  shown  for 
MLD  standard  deviations  in  figures  4.5  to  4.8. 

These  plots  show  that  most  of  the  variability,  for 
all  definitions,  is  associated  with  the  change  of  season  and 
r el  at  i  v  e  1  y  1 i 1 1 1  e  is  associated  with  a  change  of  area. 
On  1  y  for  mLDjh  do  these  statistics  show  that  the  mean 
of  area  II  is  deeper  than  the  mean  of  area  I.  The  other 
def i ni iions  show  only  that  the  mean  of  area  II  is  deeper 
than  the  mean  of  area  III.  Th at  TT~,  the  95 i  confidence 
intervals  for  these  cases  do  not  overlap. 

Only  for  the  density  and  stability  criteria  is 
the  deviation  Tn  area  IT  significantly  h i gher  than  the 
dev  i  at  i  on  i  n  area  I.  These  plots  do  show  that  in  fall 
the  deviations  in  Area  III  are  larger  than  the  deviations 
in  Area  I  for  aTl  criter i a  except  stabi 1 ity.  Tn  general , 
more  appropriate  statistical  tests  are  required  to  demon¬ 
strate  an  area  dependence  in  MLD  distributions.  Non- 


parametrlc  statistics  will  provide  the  necessary  tests 
in  the  following  section. 

The  rather  high  uncertainty  In  spring  and  fall 
for  all  definitions  Is  caused  by  the  large  standard  devi¬ 
ations  of  those  distributions  which  themselves  are  caused  by 


TABLE  4.1 

MEANS  AND  STANDARD  DEVIATIONS  OF  MLD  DISTRIBUTIONS 

(in  meters) 


95*  95*  95*  95* 
CLASS--  MEAN  LOWER  UPPER  STANDARD  LOWER  UPPEI 
CRITERION  MLD  BOUND  BOUND  DEVIATION  BOUND  BOUNI 


Summer,  Area  I 


T 

r 

V 

D 

N 

28.5 
21.8 
35.0 

15.5 

21.9 

14.9 
28.0 
12.6 

35.1 
28.7 

42.1 
18.4 

22.2 

23.2 

22.3 

9.8 

18.4 
19.3 

18.5 
8.1 

28.0 

29.2 

28.3 

12.3 

Fall,  Area 
T 

I 

65.1 

55.5 

74.8 

12.6 

8.5 

24.0 

C 

57.1 

43.4 

70.9 

17.9 

12.1 

34.3 

D 

73.3 

65.6 

81.0 

10.0 

6.8 

19.2 

N 

28.7 

10.5 

46.8 

23.6 

15.9 

45.2 

Summer,  Area  II 

T  36.6 

32.7 

40.5 

17.3 

14.9 

20.5 

C 

30.3 

26.5 

34.1 

16.9 

14.6 

20.0 

D 

63.7 

32.7 

94.7 

136.4 

117.8 

162  1 

N 

22.9 

19.4 

26.5 

15.7 

13.6 

18.7 

Fall,  Area 
T 

ir 

76.2 

51.7 

100.8 

32.0 

21.6 

61.2 

C 

69.4 

44.9 

94.0 

31.9 

21.6 

61.2 

D 

85.9 

66.4 

105.3 

25.3 

17.1 

48.5 

N 

47.2 

24.4 

70.1 

29.7 

20.1 

56.9 

Spring,  Area  III 

T  53.6 

26.9 

80.3 

74.0 

59.4 

98.4 

C 

65.9 

15.3 

116.6 

140.5 

112.6 

186.7 

D 

138.5 

72.7 

204.4 

179.5 

141.6 

234.7 

N 

27.6 

13.2 

42.0 

39.9 

32.0 

53.1 

Summer,  Area  III 

T  23.0 

20.9 

25.0 

8.3 

7.1 

10.0 

C 

16.5 

14.5 

18.6 

8.2 

7.0 

9.9 

0 

27.6 

24.7 

30.6 

11.7 

9.9 

14.2 

N 

18.6 

15.9 

21.4 

10.9 

9.3 

13.2 
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TABLE  4.1  (continued) 

MEANS  ANO  STANOARD  DEVIATIONS  OF  MLD  DISTRIBUTIONS 

(in  meters) 


CLASS-- 

MEAN 

95% 

LOWER 

95% 

UPPER 

STANDARD 

95% 

LOWER 

95% 

UPPER 

CRITERION 

MLD 

BOUND 

BOUND 

DEVIATION 

BOUND 

BOUND 

Fall,  Area  III 


T 

C 

D 

N 

75.8 

132.2 

162.7 

33.6 

45.1 
79.5 

48.2 
12.4 

106.6 

184.9 

277.2 

54.9 

45.8 

78.5 
160.1 

31.6 

32.0 

54.8 

110.1 

22.1 

80.3 
137.7 
292.3 

55.4 

All 

Area  I 

T 

34.5 

27.8 

41.3 

24.9 

21.0 

30.7 

C 

27.6 

20.6 

34.6 

25.9 

21.8 

31.9 

D 

41.5 

34.6 

48.5 

25.3 

21.2 

31.3 

N 

17.6 

14.0 

21.3 

13.7 

11.5 

16.8 

All 

Area  II 

T 

40.8 

36.0 

45.5 

22.6 

19.7 

26.6 

C 

34.4 

29.7 

39.1 

22.3 

19.4 

26.19 

0 

66.1 

38.7 

93.4 

129.4 

112.5 

152.24 

N 

25.5 

21.4 

29.5 

19.0 

16.5 

22.3 

All 

Area  III 

T 

37.7 

28.7 

46.7 

47.1 

41.5 

54.3 

C 

43.4 

26.6 

60.3 

88.7 

78.2 

102.3 

D 

74.1 

50.3 

97.9 

123.1 

108.3 

142.7 

N 

22.9 

18.0 

27.8 

25.8 

22.8 

29.7 

All 

Stations 

T 

38.1 

33.6 

42.5 

35.6 

32.8 

39.2 

C 

36.8 

29.2 

44.4 

61.0 

56.2 

67.2 

D 

64.1 

50.0 

78.3 

112.3 

103.5 

123.9 

N 

22.6 

20.0 

25.3 

21.4 

19.8 

23.6 
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Figure  4.1  Dependence  of  MLD  on  area  and  season  for 
temperature  criterion.  Mean  MLDT  and  95% 
confidence  interval  about  mean.  Areas  are 
denoted  by  symbols  and  seasons  arranged 
along  horizontal  axis.  Areas  within  a  sea¬ 
son  are  offset  for  clarity. 
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4.3  Dependence  of  MID  on  area  and  season  fo 
density  criterion.  Wean  flL&jj  and  95*  con 
fldence  Interval. 
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Figure  4.4  Deoendence  of  MLO  on  area  and  season  for 
stability  criterion.  Mean  KLDN  and  95*  con¬ 
fidence  Interval  about  mean. 
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Dependence  of  ML D  standard  deviation  on 
area  and  season  for  the  sound  speed  criterion. 
MLD^  standard  deviation  and  95%  confidence 
Interval  on  the  standard  deviation. 
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Figure  4.7 


Dependence  of  MLD  standard  deviation  and  95% 
confidence  interval  on  area  and  season  for 
the  density  criterion.  Note  break  in  vertical 
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Figure  4.8 
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Dependence  of  MLD  standard  deviation  on 
area  and  season  for  the  stability  criterion 
MLD  standard  deviation  and  S5X  confidence 
interval  on  the  standard  deviation. 


the  occasional  recording  of  extremely  deep  MLDs  (in  excess 
of  several  hundred  meters)  combined  with  small  sample  sizes 
at  those  times  of  years. 

The  various  definitions  of  Ml D  show  the  same 
increase  in  both  mean  and  standard  deviation  in  the  non¬ 
summer  months .  Only  the  density  definition  57  RTTD  Tn  the 
summer  in  area  II  breaks  the  pattern.  The  increased  mean 
and  standard  deviation  for  this  definition  is  caused  by  two 
stations  with  large  values  for  MLDp  (one  with  a  value  of 
1070m).  Further  work  should  investigate  these  stations  for 
possible  problems  with  the  data  or  the  implementation  of  the 
density  criterion.  Although  it  must  now  be  considered 
unlikely,  the.  possibility  of  layers  actually  homogeneous  in 
density  with  such  fantastic  thickness  must  also  be  reviewed. 


4.2  GAUSSIANITY  OF  MLD  DISTRIBUTIONS 

The  chi-square  goodness-of-fit  test  is  applied 
to  determine  the  probability  that  a  particular  sample  dis¬ 
tribution  is  drawn  from  a  Gaussian  population.  Table  4.2 
presents  the  results  of  these  tests.  The  null  hypothesis, 
that  an  MLD  distribution  is  Gaussian,  will  be  rejected 
whenever  the  probability  of  occurrence  of  a  chi-square  value 
greater  than  or  equal  to  the  computed  test  statistic  is  less 
than  0.05.  Table  4.2  shows  the  significance  level  attained 
in  each  test  of  the  null  hypothesis.  Using  this  criterion 
the  distributions  of  the  following  classes  and  definitions 
can  be  considered  non-Gaussi an:  MLDy  and  ML  Dc  for  area  I, 
summer;  MLDq,  MLDq,  MLOr  for  area  II,  summer;  all  defini¬ 
tions  for  area  III,  spring;  MLD^  for  area  III,  summer;  all 
definitions  for  area  I,  all  seasons  combined;  all  defini¬ 
tions  for  area  II,  all  seasons  combined;  all  definitions  for 
area  III,  all  seasons  combined;  and  all  definitions  for  all 
area,  all  seasons  combined.  The  fall  season  distributions 
do  not  contain  enough  values  to  perform  the  chi  squared 
test.  The  d i str i but i on s  t h at  can  still  be  considered 
Gaussian  are:  ML  Dq  and  ML  tty  for  area  I,  summeT;  MlDj  for 
area  II,  summer;  and  MLDy,  MLDq  and  MLDq  for  area  III, 
summer 


4.3  MLD  DEPENDENCE  ON  SEASON,  AREA  AND  DEFINITION 

Non-Gaussian  distributions  can  be  tested  for  sig¬ 
nificant  differences  in  location  (mean  and  median)  by  use  of 
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TA8LE  4.2 

CHI  SQUARE  TEST  OF  GAUSSIANITY 


*  distribution  considered  Gaussian 


CLASS/DEFINITION 

CHI-SQUARE 

DEGREES  OF 

SIGNIFICANCE 

STATISTIC 

FREEDOM 

LEVEL 

Summer,  Area  I 

T 

30.69 

2 

<.005 

C 

59.20 

2 

<.005 

0 

0.64 

2 

.82  * 

N 

3.70 

2 

.12  * 

Fall,  Area  I 

i nsuf f i c i ent 

data 

Summer,  Area  II 

T 

7.23 

4 

>.10  * 

C 

10.43 

4 

.03 

D 

158.94 

4 

<.005 

N 

16.34 

4 

<.005 

Fall,  Area  II 

insufficient 

data 

Spring,  Area  III 

T 

47.60 

2 

<.005 

C 

61.70 

2 

<.005 

D 

44.19 

2 

<.005 

N 

35.94 

2 

<.005 

Summer,  Area  III 

T 

3.12 

4 

.65  * 

C 

5.17 

4 

.41  * 

D 

3.34 

4 

.63  * 

N 

10.12 

4 

.03 

Fall,  Area  III 

i nsuf f ic i ent 

data 
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TABLE  4.2  (continued) 


CLASS/OEFI N IT  ION 

CHI  SQUARE 
STATISTIC 

DEGREES  OF 
FREEDOM 

SIGNIFICANCE 

LEVEL 

All  Seasons, 

Area  I 

T 

6.69 

2 

.01 

C 

13.40 

2 

<.005 

D 

6.58 

2 

.01 

N 

6.80 

2 

.01 

All  Seasons, 
Area  II 


T 

13.12 

4 

.01 

C 

18.79 

4 

<.005 

D 

126.69 

4 

<.005 

N 

20.09 

4 

<.005 

All  Seasons, 
Area  III 


T 

138.23 

4 

<.005 

C 

166.42 

4 

<.005 

D 

229.30 

4 

<.005 

N 

95.30 

4 

<•005 

All  Seasons, 

All  Areas 

T 

175.63 

4 

<.005 

C 

269.27 

4 

<.005 

D 

305.36 

4 

<.005 

N 

83.78 

4 

<.005 

1 


non-parametrlc  statistical  methods.  These  methods  do  not 
depend  upon  the  characteristics  of  the  probability  density 
functions  of  the  distributions.  The  differences  in  MLD 
distributions  from  different  areas  and  seasons  can  be  tested 
using  the  Mann-Whitney  two  sample  test,  which  does  not 
require  an  equal  number  of  observations  in  each  sample.  MLD 
distributions  based  on  different  definitions  were  compared 
using  the  Wilcoxon  Signed  Rank  Test.  The  Hodges-Lehmann 
estimator  was  used  for  differences  in  the  mean  MLD  and  the 
associated  confidence  interval  for  the  difference  in  the 
mean  MLD. 


The  seasonal  dependence  of  MLD  is  examined  first. 
The  Mann  Whitney  test  results  are  given  in  table  4.3.  The 
significance  level  indicates  the  probability  that  the  test 
statistic  could  be  produced  from  two  distributions  which 
come  from  the  same  population.  When  this  value  drops  to  5% 
or  less,  then  the  distributions  are  considered  significantly 
different.  Those  d istribut  ions  not  considered  different  are 
marked  by  a  double  asterisk,  while  nearly  different  cases 
are  marked  with  a  single  asterisk.  The  table  also  contains 
an  estimate  of  the  difference  in  the  means,  alonq  with  upper 
and  lower  95<  confidence  limits  on  the  estimate.  Plots  of 
the  estimated  difference  in  the  means  and  the  95X  confidence 
intervals  are  presented  in  figure  4.9a,  b,  and  c. 

Summer  and  fall  distributions  (figure  4.9a)  are 
significantly  diYTerent  for  all  areas  "and  for  all  criterTa 
except  N.  The  difference  is  about  45  +  10  m  in  areas  I  and 
II  (fall  being  deeper)  where  it  is  independent  of  criterion 
(excluding  N).  The  difference  is  greater  in  area  III  (about 
70  +  29  m),  more  variable  and  more  dependent  on  criterion. 
In  summer,  area  III  has  the  shallowest  mean  MLD;  in  winter, 
the  deepest  (see  figures  4.1  to  4.3).  How  this  characteris¬ 
tic  may  be  related  to  frontal  zone  dynamics  in  the  Iceland 
Faroe  Front  is  not  clear. 

The  lack  of  a  seasonal  variation  in  MLDn  distri¬ 
butions  points  to  a  serious  prop i em  with  that  criterion.  It 
Ts  likely  that  noise  still  in  the  data  is  generating  layers 
of  high  N  values  where  none  actually  exist  (see  the  struc¬ 
ture  still  apparent  in  the  smoothed  o j  profile,  figure 
2.5).  This  suqgests  that  smoothing  to  only  5  m  is  not  suf¬ 
ficient  to  remove  all  the  time  lag  problems  in  the  data. 
Such  noise  will  tend  to  trlqger  the  MLD#  algorithm  prema¬ 
turely,  and  hence  explain  the  shallow  nature  and  seasonal 
independence  of  the  MLDr  distributions. 
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TA8LE  4.3 

SEASONAL  DIFFERENCES  BY  MANN  WHITNEY  TESTS 
*  marginally  insignificant  difference 
**  insignificant  difference 


SEASON,  AREA 
CRITERION 

SIGNIFICANCE 

LEVEL 

MEAN 

DIFFERENCE 

95% 

LOWER 

BOUND 

95% 

UPPER 

BOUND 

Fal 1-Summer, 

T 

Area 

I 

<.0014 

41 

30 

49 

C 

<.0014 

42 

28 

51 

0 

<.0014 

42 

26 

53 

N 

.1038** 

7 

-2 

30 

Fal 1 -Summer , 

T 

Area 

II 

<.0014 

50 

28 

62 

C 

<.0014 

49 

29 

62 

D 

<.0014 

48 

35 

59 

N 

.0038 

16 

5 

38 

Fal 1-Summer, 

T 

Area 

III 

<.0014 

67 

13 

88 

C 

<.0014 

93 

79 

115 

D 

<.0014 

74 

33 

111 

N 

.2843** 

6 

-6 

27 

Spr 1 ng-Summer 
T 

,  Area  III 

<.0014 

9 

3 

16 

C 

.0055 

7 

1 

14 

0 

<.0014 

19 

7 

99 

N 

.1230** 

2 

-2 

7 

Fal  1-Spr ing, 

T 

Area 

III 

.0475* 

45 

-3 

75 

C 

<.0014 

84 

60 

103 

0 

.2148** 

28 

-42 

94 

N 

.4443** 

2 

-9 

25 
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The  summer-spring  comoarlson  (figure  4.9b)  shows 
significant  differences  for  temperature  and  density,  and  a 
nearly  significant  difference  for  sound  speed.  As  before, 
stability  does  not  respond  to  the  seasonal  change.  Spring 
MLOs  are  about  10m  deeper,  on  average,  than  summer  MLDs. 

The  spring-fall  differences  are  significant  for 
the  sound  speed  criterion  and  marginally  significant  for  the 
temperature  criterion.  The  density  distribution  does  not 
show  a  significant  difference.  As  before,  the  MLDr  values 
are  insensitive  to  seasonal  differences. 

The  dependence  of  MLD  on  area  is  examined  next. 
The  Mann  Whitney  test  results  are  given  in  table  4.4,  and 
the  mean  differences  and  95%  confidence  limits  are  presented 
in  figures  4.10a,  b  and  c. 

In  summer,  area  II  is  significantly  deeper  than 
area  I  for  all  criteria  (figure  4.10a).  Large  confidence 
intervals  are  associated  with  the  mean  difference  of  the 
density  MLDs  because  there  is  a  great  deal  of  variability 
in  the  differences  due  to  the  deep  observations  of  MLD 
in  area  II  in  the  summer  under  that  criterion.  In  fall, 
area  II  is  significantly  deeper  than  area  I  for  temperature 
and  stability  and  nearly  significantly  for  sound  speed. 
Area  II  is  deeper  than  area  I  by  about  10m  in  the  summer 
aind  20m  in  the  fall. 

Area  III  distributions  of  MLD  cannot  be  considered 
different  from  those  of  area  I  for  any  criterion  in  summer 
or  any  criterion  but  sound  speed  in  fall. 

In  summer,  area  III  is  significantly  shallower 
than  area  II  for  all  criteria  by  an  approximate  value  of 
10m.  In  fall,  only  the  sound  speed  gives  a  significant 
difference.  In  this  case,  area  III  is  deeper  than  area  II. 

Finally,  the  dependence  of  MLD  on  criterion  is 
examined.  The  temperature  definition  is  taken  arbitrarily 
as  a  standard. 

The  differences  between  the  MLD  means  for  the 
sound  speed  and  temperature  criteria  as  a  function  of  season 
and  area  are  presented  In  figure  4.11a.  Except  for  fall  in 
area  III,  the  MLDq  values  tend  to  be  6  m  +  2  m  shallower 
than  MLDj  values.  The  discrepancy  in  area  III  during  the 
fall  Is  significant  and  would  require  more  study  to  relate 
to  oceanographic  properties. 
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TABLE  4.4 

AREA  DIFFERENCES  BY  MANN  WHITNEY  TESTS 
*  marginally  Ins ign 1 f i cant  differences 
**  insignif icant  differences 


95* 

95* 

AREA, 

SEASON 

SIGNIFICANCE 

MEAN 

LOWER 

UPPER 

CRITERION 

LEVEL 

DIFFERENCE 

BOUND 

BOUND 

Area 

II-Area  I 

,  Summer 

T 

<.0014 

11 

5 

16 

C 

<.0014 

11 

5 

17 

D 

.0080 

10 

-3 

23 

N 

.0019 

6 

1 

10 

Area 

I  I-Area  I 

,  Fall 

T 

.0351 

22 

-26 

37 

C 

.0516* 

22 

-11 

36 

* 

0 

.0122 

19 

6 

32 

N 

.0721** 

15 

-17 

49 

Area 

Ill-Area 

I,  Summer 

T 

.3557** 

0 

-5 

4 

C 

.5120** 

0 

-3 

4 

0 

.0594* 

-5 

-12 

2 

N 

.0853** 

3 

0 

7 

Area 

Ill-Area 

I,  Fall 

T 

.1922** 

26 

-37 

53 

C 

<.0014 

52 

27 

92 

D 

.1093** 

23 

-39 

267 

N 

.4840** 

2 

-21 

35 

Area 

I  II-Area 

II,  Summer 

T 

<.0014 

-13 

-17 

-8 

C 

<.0014 

-12 

-17 

-7 

D 

<.0014 

-20 

-30 

-1  1 

N 

.0436 

-3 

-7 

1 

Area 

Ill-Area 

II,  Fall 

T 

.3121** 

7 

-55 

32 

C 

.0044 

37 

11 

85 

D 

.3121** 

15 

-53 

249 

N 

.1190** 

-16 

-35 

13 
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The  differences  in  MLD  mean  values  for  the  density 
and  temperature  criteria  as  a  function  of  area  and  season 
are  presented  in  figure  4.11b.  For  areas  I  and  II,  the 
MLOq  means  are  about  5  m  deeper  than  the  MLD7  means. 
For  area  III,  the  difference  in  the  means  is  much  greater 
and  much  less  certain.  Stated  another  way.  In  area  III  the 
density  criterion  MLDs  respond  much  more  to  the  sea¬ 
sonal  cycle  than  the  temperature  MLDs. 

The  difference  between  temperature  and  stability 
(figure  4.11c)  shows  once  again  the  failure  of  the  MLDN 
means  to  respond  to  the  seasonal  cycle. 

Overall  (figure  4.12),  the  sound  speed  MLDs  do 
not  differ  from  the  temperature  MLDs,  the  density  MLDs  are 
deeper  and  the  stability  MLDs  are  shallower. 
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Figure  4.9 


Seasonal  dependence  of  MLD  as  a  function  of 
area  and  criterion.  Symbols  represent  cri¬ 
teria  and  position  on  horizontal  axis  repre¬ 
sents  area.  Within  an  area  symbols  are  offset 
for  clarity.  a)  fall-summer. 


Figure  4.10  Area  dependence  of  MLD  as  a  function  of 
season  and  criterion.  Symbols  denote 
criteria  and  horizontal  position  repre¬ 
sent  season.  Within  a  season  symbols  are 
offset  for  clarity,  a)  Area  II  -  Area  I 
b)  Area  III  -  Area  I. 


F/G  20/1 


AD-A105  669 
UNCLASSIFIED 


SCIENCE  APPLICATIONS  INC  MCLEAN  VA 

MIXED  LAYER  DEPTH  DISTRIBUTIONS  FROM  SVSTD  DATA. IU) 

APR  81  E  MOL I NELL I »  J  OONELSON.  L  LILLY  N6B463-80 

SAI-81-260-WA  N00-TR-271 


-C-0035 

NL 


SPRING 


SUMMER 


gure  4.11 


a  Difference  between  sound  speed  and 
temperature  criteria  as  a  function 
of  area  and  season.  Mean  difference 
and  95%  confidence  interval. 


SPUING 


SUMMER 


FALL 


■  ALL  AREAS 
ALL  SEASON; 


Figure  4.12  Dependence  of  MLD  on  criterion.  All  data  for 
all  seasons  and  areas  combined.  The  sound 
speed  criterion  Is  not  significantly  different 
than  the  temperature.  At  the  .001  level  of 
significance  the  MLDqs  are  deeper,  and  the 
MLDrs  are  shallower,  than  the  MLDts. 


Section  5 
CONCLUSIONS 


Five  conclusions  can  be  drawn  from  this  work. 

1.  The  greatest  variation  of  MLD  occurs  with 

changes  of  season  and  amounts  to  an  Increase 
of  40m  or  more  from  summer  to  fall. 

2.  There  is  relatively  slight  but  still  signifi¬ 

cant  variation  of  MLD  with  area.  The  Iceland 
Flank  MLOs  are  about  10m  deeper  than  either 
the  Rockall  Trough  or  Iceland  Faeroe  Frontal 
Zone  in  the  summer  and  perhaps  even  more  In 
the  fall.  There  Is  not  a  significant  differ¬ 
ence  between  the  Rockall  Trough  and  the 

Iceland  Faeroe  Frontal  Zone. 

3.  The  definition  of  MLD  based  on  stability  is 
not  sensitive  to  the  seasonal  variations,  but 
seems  to  be  affected  by  noise  in  the  density 
profile  introduced  by  time  lag  problems  and 
not  removed  by  5m  filtering. 

4.  Except  for  Area  III  in  the  spring  and  fall, 
the  difference  between  the  other  definitons  of 
MLD  are  relatively  independent  of  area  and 
season.  Area  III  is  also  remarkable  in  having 
more  deviations  in  MLD  than  Area  I  in  the 
fall . 

5.  The  large  values  of  MLD  in  summer  in  Area  II 
are  suspect.  These  occur  at  stations  15  and 
18  of  cruise  932014. 


Section  6 
RECOMMENDATIONS 


Based  on  the  work  performed  here  three  tasks  for 
further  work  can  be  identified.  The  first  Improves  the 
present  analysis,  the  second  attempts  to  generalize  the 
results  and  the  third  applies  to  future  survey  planning. 


6.1  IMPROVEMENTS 

Some  of  the  MLD  criteria  should  be  redefined 
to  be  more  consistent  with  applications. 

The  present  definition  of  the  stability  MLD 
should  not  be  used  in  future  analyses  as  it  is  too  sensitive 
to  noise  in  the  data.  All  MLD#  distributions  should  be 
recalculated  with  a  new  definition  that  is  triggered  by 
Brunt-Val sal  a  peaks  in  excess  of  3  cph  (2  cph  is  presently 
used) . 


The  definition  of  the  sound  speed  MLD  should 
take  into  account  the  sound  speed  gradient  above  the 
maximum.  Without  a  sufficient  Increase  of  sound  speed  with 
depth,  surface  ducting  of  sound  energy  will  not  occur.  A 
sensitivity  test  of  sound  energy  ducting  as  a  function  of 
sound  speed  gradient  in  the  mixed  layer  must  be  performed. 
Then  MLD^  should  be  recalculated  with  a  requirement  that 
gradient  above  the  sound  speed  maximum  exceed  a  minimum 
value.  Then  the  resulting  distribution  of  MLDc  will 
pertain  to  surface  ducting. 

The  present  analysis  will  also  be  Improved 
by  simply  eliminating  bad  stations.  Two  summertime  stations 
have  been  Identified  as  suspect.  These  should  be  examined 
more  closely  to  look  for  causes  for  the  discrepancy  between 
their  MLD  values  and  the  values  of  MLD  from  other  stations 
In  the  area.  If  unrealistic  measurements  are  to  blame  the 
stations  should  be  dropped. 

In  order  to  better  characterize  a  particular 
distribution  of  MLDs  (only  MLDf  need  be  characterized  If 
the  relationships  with  the  other  criteria  are  known,  see 
below)  more  needs  to  be  known  than  the  fact  that  the  dls- 
trlbutlon  Is  non-GaussI an .  By  plotting  the  third  and  fourth 
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moments  of  the  sample  distribution,  an  analytical  distri¬ 
bution  (Binomial,  Poissan,  etc.)  can  be  selected  that  Is 
most  appropriate  for  the  data.  This  last  step  does  not  lead 
to  a  more  complete  understanding  of  the  processes  leading  to 
the  distributions  of  MLDs,  it  just  makes  a  more  complete 
description  easier  to  specify. 


6.2  GENERALIZING  RESULTS 

It  would  be  very  useful  if  MLDs  by  the  various 
definitions  could  be  predicted  from  the  values  of  MLD 
according  to  the  temperature  criterion.  Temperature  profiles 
are  both  simpler  and  much  less  expensive  to  obtain  (an 
expendable  bathythermograph,  XBT,  can  be  used  instead  of 
an  SVSTD).  XBT  measurements  are  more  plentiful  in  the 
historical  data  set.  For  this  reason  It  is  highly  recom¬ 
mended  that  more  research  be  performed  to  correlate  MLDf 
to  MLDc,  MLDq  and  MLDs*  The  data  set  provided  for  the 
present  study  is  already  useful  to  make  the  correlations  for 
summer  conditions. 

The  relationship  can  be  investigated  using 
three  approaches.  First,  look  at  the  equations  of  state 
relating  sound  speed,  density  and  density  gradient  to 
temperature  (T)  and  salinity  (S)  to  determine  how  the 
various  criteria  should  be  related  for  a  given  T-S  rela¬ 
tionship.  Assume  some  simple  T-S  curves  and  predict  the 
correlation  o7  MlDj  to  MLDc,  ML  Do  and  MLDs .  Second , 
use  the  S  T  D  d  at  a  Tn  the  area  to  define  the  actual  T-S 
curve(s)  and  use  these  to  calculate  a  second  set  of  corre- 
1  ations  between  mU5j  and  MLDc,  ML()d  and  MLDs-  See  Tf 
these  differ  significantly  from  the  correlations  computed 
using  the  simplified  T-S  relations.  Finally,  compute  the 
actual  correlations  between  MLDj  and  MLftc ,  MLOq  and  ML  Os 
and  see  T7  these  are  consistent  with  the  eguations  of 
state  and  the  observed  or  simplified  T-S  curves. 


6.3  FUTURE  SURVEY  PLANS 

To  make  useful  correlations  in  other  seasons 
more  data  must  be  taken.  There  are  statistical  methods  to 
determine  the  number  of  stations  that  must  be  occupied  in  a 
particular  time-space  window  in  order  to  characterize  the 
distribution  of  some  property  (e.g.  MLD)  in  that  window. 


96 


In  addition,  statistical  methods  exist  that 
can  provide  sampling  strategies  to  address  area  charac¬ 
terizations.  That  Is,  they  address  the  question  of  how 
large  the  time-space  window  mentioned  above  should  be. 
These  methods  typically  rely  on  estimates  of  correlation 
lengths  and  correlation  times.  Some  of  these  correlation 
scales  can  be  determined  from  individual  cruises  among  the 
data  provided  for  the  present  study.  Correlation  times, 
specifically,  will  require  additional  data  sets,  preferrably 
from  moored  Instruments  as  used  in  the  JASIN  area. 


Appendix  A 

COMPUTER  PROGRAMS  FOR  MLD  DEFINITIONS 


This  appendix  lists  the  FORTRAN  statements  that 
are  used  to  select  MLDs  from  arrays  of  pressure  (P,  in  dbar) 
temperature  (T,  in  °C)  and  salinity  (S,  in  0/00  for  the 
calculation  of  sound  speed,  density  and  Brunt- V*ii  sill  fre¬ 
quency).  In  these  routines,  N  is  the  number  of  elements  in 
the  P,  T  and  S  arrays.  The  value  of  mixed  layer  depth  is 
returned  (TMLD ,  SSMLD ,  OMLD  or  BVMLD)  as  is  the  index  of  the 
array  element  (IMLD)  associated  with  that  depth.  For  the 
calculation  of  BVMLD,  a  value  of  the  local  gravitational 
acceleration  is  passed  as  the  variable  G.  Jhe  calculations 
of  sound  speed,  density  and  Brunt-Vai sal  a  frequency  are 
detailed  in  section  2  of  this  report. 
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A .  1 


FORTRAN  SUBROUTINE  FOR  MLDj 


SUBROUTINE  TEMMLO  { P ,  T,  S,  N,  TMLD,  IMLD) 
DIMENSION  P(l),  T(  1 ) ,  S(l) 

If  ( (P(l)-5.0) .GE.0.5)  GO  TO  500 
DO  10  1*1, N 
C  FIND  5  M  FIRST 

IF  (P(I)-5. .GE.0.5)  GO  TO  20 

15*1 

T5*T( I ) 

GO  TO  10 

C  CHECK  AGAINST  T  VALUE  AT  5  DB 
20  TDI F*ABS( T( I ) -T5 ) 

IF  (TDIF.LT.0.2)  GO  TO  10 
TGRAD*ABS( (T( I)-T(I-l) )/(P( I)-P(I-l) )) 

C  IF  TGRAD  LE  5.E-3  KEEP  LOOKING 

IF  (TGRAD. LE . 5 . E-3 )  GO  TO  10 
C  STOP----DESIRED  VALUES  FOUND 

TMLD*P( I) 

IMLD* I 
RETURN 
10  CONTINUE 
500  TMLD*-10 . 

IML0-00 

RETURN 

END 
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.2 


FORTRAN  SUBROUTINE  FOR  ML Dc 


SUBROUTINE  SMLD(  P , T , S , N , SSMLD, IMLD) 

OIMENSION  P(  1 ) ,  T(  1 )  ,  S(  1 ) 

DIMENSION  SS  (20000) 

ERR* .  5 

IF((P(l)-5.0) .GE.0.5)  GO  TO  500 
M  =  1 

SSMAX *0 . 0 
00  10  1*1, N 

I F ( P ( I ) .LT.5.0)  GO  TO  10 

ROUND  SOUND  SPEEO  BECAUSE  ADJ  VALUES  MAY  BE  SIMILAR 

IS*(SOUND(P(  I) ,T(I) ,S(I) )+ . 05 ) *10 . 0 

SS( I ) *  I  S/10 .0 

CHECK  FOR  A  NEW  MAXIMUM 

IF(SS(I) .LT. SSMAX)  GO  TO  95 

SSMAX*SS( I ) 

M=  I 

GO  TO  10 
95  CONTINUE 

SSDIFF=SSMAX-SS(I) 

CHECK  WHETHER  VALUES  FOUND 
I F( SSD IFF .GE . ERR)  GO  TO  99 
GO  TO  10 
99  SSMLD*P( M) 

IMLO=M 
RETURN 
10  CONTINUE 
500  CONTINUE 
SSMLD* -10 . 

IMLD*0 

RETURN 

END 


FORTRAN  SUBROUTINE  FOR  MLO, 


A. 3 


o  T 


SUBROUTINE  DENMLO(P,T,S,N,DMLO,IMLD) 

DIMENSION  P(l) ,T(1) ,S(1) 

IF( (P(l)-5.0) . GE . 0 . 5 )  GO  TO  500 
15-1000 
DGMAX-0 .0 
DO  101*1  ,N 

IF( (P(I)-5.0) .GE.  .5)  GO  TO  20 
15*1 

CALL  SIGGT{T(I),S(I),S0,ST5) 

ST0LD-ST5 
GO  TO  10 

20  CALL  SIGGT(T( I ) , S( I ) , SO , ST) 

C  CHECK  AGAINST  SIGMA  T  VALUE  AT  5  DBARS  TO  SEE  IF  FOUND 
IF( (ST-ST5) .LT.0.02)  GO  TO  99 
OG-ST-STOLD 

IF(OG.LE.OGMAX)  GO  TO  30 

DGMAX-DG 

IDGMX-I 

GO  TO  99  - 

30  IF(DGMAX.LE.O.O)  GO  TO  99 
DGDIF-DGMAX-DG 

C  IF  OGOIF  LT  ERROR  IN  DG  THEN  KEEP  LOOKING 
IF(DGDIF.LE.O.Ol)  GO  TO  99 
IMLO-IOGMX 
DMLD-P(IMLD) 

RETURN 
99  STOLD-ST 
10  CONTINUE 
500  0ML0--10 . 

IMLO-O 

RETURN 

END 
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A. 4  FORTRAN  SUBROUTINE  FOR  MLDfj 


SUBROUTINE  BMLD(P,T,S,N,BVMLD, IMLD.G) 
DIMENSION  P(1),T(1),S(1) 

DIMENSION  B V ( 2000) , PNSQ( 2000) 

I F( N . GT. 2000)  N*2000 
ERR* .  5 

IF(P(l)-5.  .6E. .5)  SO  TO  500 
PNSQ( 1)*2 .0 

CALL  BVSQ(P,T,S,S,N,BV,PNSQ) 

BVMAX*0 . 

I BMX*0 
DO  10  I  *1 ,  N 

I F ( B V ( I) .LT.O.O)  GO  TO  199 
BV ( I ) *SORT( BV( I )  ) 

GO  TO  10 
199  B V( I ) *0 .0 
10  CONTINUE 
00  200  1*1, N 
C  FIND  BV  MAX  VALUE 

I F ( B V ( I) .LE.2.0.AND.BVMX.LE.0.0)  GO  TO  200 
IF( P( I ) . LT. 5 .0 )  GO  TO  200 
IF(BV( I) .LE.BVMAX)  GO  TO  30 
BVMAX*BV( I) 

I8MX*I 
GO  TO  200 

30  IF((BVMAX-BV( I) ) .LE.ERR)  GO  TO  200 
IMLD* I BMX 
BVMLD*P( IMLD) 

RETURN 

200  CONTINUE 
500  BVMLD*-10. 

IML0=0 

RETURN 

END 
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ACRONYMS 


ML  D 

NAVOCEANO 

SAI 

SVSTD 


mixed  layer  depth 

Naval  Oceanographic  Office 

Science  Applications,  Inc. 

sound  velocity,  salinity, 
temperature,  depth  In  situ 
continuously  recording  probe 
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SYMBOLS 


OC 

o/oo 
Area  I 
Area  II 
Area  III 
C 

Cl 

D 

9 

N 

P 

PP 

PR 

S 

T 

Z 

a 

Y 

e 

p 

♦ 


degrees  Celsius 
parts  per  thousand 
Rockall  Trouqh 
Iceland  Flank 

Iceland  Faroe  Frontal  Zone 

sound  speed 

chlorlnlty 

density 

gravitational  acceleration 
Brunt-V*41  s‘411  Frequency 
pressure 

effective  pressure  for  adiabatic 

lapse  rate  calculation 

reference  pressure 

salinity 

temperature 

depth 

in  situ  specific  volume 
171 ab at  1c  lapse  rate 
potential  temperature 
density 

specific  qravlt.y  anomaly 
1 atitude 
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